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We show that the colour glass condensate dipole model of lancu, Itakura and Munier, improved 
to include the impact parameter dependence, gives a good fit to the total 7*p cross section measured 
at HERA if the anomalous dimension at the saturation scale, 7^, is treated as a free parameter. We 
find that the optimum value of 7s = 0.46 is close to the value determined from numerical solution 
of the Balitsky-Kovchegov equation. The impact parameter dependent saturation scale is generally 
less than 0.5 GeV^ in the HERA kinematic regime for the most relevant impact parameters b ~ 2-3 
GeY-\ We compare predictions of the model to data on the longitudinal and heavy flavour structure 
. . . functions, exclusive diffractive vector meson production and deeply virtual Compton scattering at 

00 ' HERA. The model is found to be deficient for observables sensitive to moderately small dipole sizes, 

, where an alternative model with explicit DGLAP evolution performs better. The energy dependence 

^— »' ' of exclusive diffractive processes is shown to provide an important discriminator between different 

dipole model cross sections. 
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I. INTRODUCTION 



The colour dipole model has proven to be very successful in describing a wide variety of small- a; inclusive and 
diffractive processes at HERA. In particular, it is commonly used in determinations of the saturation scale, that is, 
' the x-dependent momentum scale at which nonlinear effects start to become important. The total cross section for 
II ■ 7*p scattering, 

^ : ^7o? = 4' + -1" = ^2, (1) 

, is obtained by combining the light-cone wave functions for the virtual photon to fluctuate into a qq pair with the 
J> ' dipole cross section for the qq pair to scatter elastically off the proton: 

\ Here, / is the flavour of the qq pair, z is the fraction of the photon's light-cone momentum carried by the quark, 
r = \r\ is the transverse size of the qq dipole, while h is the impact parameter, that is, h is the transverse distance 
from the centre of the proton to the centre of mass of the qq dipole. The squared photon wave functions, (5'*^')^^, 
are given explicitly in Ref. 1] . 

A popular parameterisation for the 6-integrated dipole cross section is the "saturation model" due to Golec-Biernat 
and Wiisthoff (GBW) [2,1]: 

H : 

^d^6^=.„(l-e-^'?=(^)/^), (3) 

where ctq is a constant and Q\(x) — [xq/x)^ GeV^. The parameters ctq = 29 mb, A = 0.28 and xq = 4 x 10^^ were 
obtained from a fit, including charm quarks, to inclusive deep-inelastic scattering (DIS) data ^2]. However, the GBW 
model does not give a good fit to recent DIS data 

Many improvements to the seminal but simple GBW parameterisation ^ of the dipole cross section have been 
proposed (see Ref. [l| for a more comprehensive review and further references). The GBW model was modified 
to include DGLAP evolution in the Bartels-Golec-Biernat-Kowalski (BGBK) model [J] (extended to include heavy 
quarks in Ref. Q). An alternative colour glass condensate (CGC) model inspired by the Balitsky-Kovchegov (BK) 
equation^ SBlSl was proposed by lancu, Itakura and Munier [loj and extended to include charm quarks in Ref. 
The introduction of charm quarks in the CGC model led to a dramatic decrease in the saturation scale [1]. However, 



^ It has recently been shown Q that running coupling effects strongly suppress the effect of "Pomeron loops" such that the BK equation 
should be sufficient for phenomenological studies. 
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it has recently been shown by Soyez [ll!| that allowin g th e anomalous dimension at the saturation scale to increase 
from the fixed value of 7^ = 0.63 assumed in Refs. [ll llOj] to a higher value of 7s — 0.74 slightly improves the fit to 
the F2 data and gives a larger saturation scale. On the other hand, this value of 7^ seems to be inconsistent with the 
value of ~ 0.44 recently obtained from numerical solution of the BK equation by Boer, Utermann and Wessels p^ . 

The BGBK 4] and CGC 10] models considered only the dipolc cross section integrated over the impact parameter 
b. However, the gluon density is larger in the centre of the proton (b = 0) than at the typical impact parameters 
b ^ 2-3 GeV~^ probed in the total 7*p cross section. Therefore, any serious determination of the saturation scale 
at HERA should consider the impact parameter dependence of the dipole cross section. The b dependence is also 
necessary to describe the slope of t distributions of diffractive processes at HERA, which in turn fix the normalisation 
of the 6-integrated dipole cross section. Therefore, in the context of dipole models, the analysis of inclusive HERA 
data should not be considered in isolation to diffractive HERA data. 

The BGBK model was modified to include the impact parameter dependence in Refs. denoted by the "b-Sat" 

model, and was found to give a good description of both F2 data and exclusive diffractive processes. The CGC model 
was also modified to include the impact parameter dependence in Ref. [H, denoted by the "b-CGC" model, but it 
was not possible to get a good description of the F2 data with a fixed value of 7^ = 0.63. It is therefore interesting 
to see if a good fit with the b-CGC model can be obtained when 7^ is allowed to vary and to compare the optimum 
value obtained with the expected value of 0.44 obtained from numerical solution of the BK equation p^ . 

In Sec. [H] we recall the original CGC model 10] and investigate the dependence on 7^ of the quality of the fit to 
HERA F2 data. In Sec. lIIII we introduce the impact parameter dependence into the model and determine the optimum 
value of 7s. The impact parameter dependent saturation scale is discussed in Sec. IIVI Predictions are obtained for 
the longitudinal and heavy flavour structure functions in Sec. |V] and for exclusive diffractive processes in Sec. IVII 
these predictions are then confronted with HERA data. Finally, we conclude in Sec. IVIII 



II. COLOUR GLASS CONDENSATE DIPOLE MODEL 

The CGC dipole cross section of lancu, Itakura and Munier ^T3], integrated over impact parameter, can be written 

as 

a,g = Jd'b^ = 2 Jd^bM{x,r,b) = 2 J d'b T{b) M{x,r) = aoMix,r), (4) 

that is, the b dependence oi J\f{x,r,b), the imaginary part of the dipole-proton scattering amplitude, is assumed to 
factorise, so that the integration over b gives a multiplicative constant ctq determined by a fit to F2 data. The usual 
interpretation [lol ] is that the proton is assumed to be a homogeneous disk of radius Rp. Then the b dependence is 
given by a step function T{b) = Q{Rp — b), so that integration over b gives ctq = 2TrRp in (g]). In fact, T{b) need not 
necessarily be a step function. For example, it is pointed out in Ref. [13] that a factorised Gaussian b dependence 
will also lead to (|4]). However, the factorisation of the b dependence from the x dependence is not supported by the 
HERA diffractive data, where one finds a significantly nonzero effective a'p, indicating correlation between the b and 
X dependence of the dipole scattering amplitude; see Sec. IVII 

The scattering amplitude J\f{x, r) can vary between zero and one, where J\f = 1 is the unitarity ("black disc") limit. 
J\f is obtained by smoothly interpolating between two limiting types of behaviour. For small dipole sizes, r <g; 2/Qs, M 
is obtained from the saddle point approximation to the leading-order (LO) BFKL equation, followed by an expansion 
to second order around the saturation saddle point. For large dipole sizes, r ^ 2/Qs, the functional form obtained 
from solving the BK equation is used [3]. The scattering amplitude is therefore [13] 

^^^ fro + ^ 

N{x,r) = \-^o[^) : rQ.<2^ 

\^l_^-Aln^{BrQ^) . > 2 

where Qs = Qs{x) = {xq/x)^/'^ GeV, Y — ln(l/a:), and k = x" {l s) / x' {l s) where x is the LO BFKL characteristic 
function. The coefficients A and B in the second line of ([5]) are determined uniquely from the condition that N{x, ?■), 
and its derivative with respect to rQs, are continuous at rQg — 2: 

For tQs ^ 2, the scattering amplitude is a function only of rQs\ this is the so-called geometric scaling. For small 
dipole sizes, rQs ^ 2, the second ("diffusion") term in the exponent of the first line of (O enhances the effective 
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Is 


(70 /mb 


Xo 


A 


xVd.o.f. 


p- value 


0.63 (fixed, [IJ) 


35.7 


2.70 X 10"^ 


0.177 


116.8/130 = 0.90 


0.79 


0.74 (fitted) 


27.4 


1.63 X 10"^ 


0.216 


110.4/129 = 0.86 


0.88 


0.61 (fitted) 


37.4 


1.09 X 10"^ 


0.170 


115.4/129 = 0.89 


0.80 


0.44 (fixed) 


46.3 


2.21 X 10'" 


0.122 


180.1/130 = 1.39 


2.4 X 10"^ 



TABLE I: Parameters of the CGC dipole model JS)), for different values of 7^, determined from fits to ZEUS F2 data [H, [13] 
with XBj < 0.01 and £ [0.25, 45] GeV'^. The meaning of the p-value is explained in the text. 



anomalous dimension^ from 7^ to 



'^°«=ain(r2Q2/4)=7s + ^ln;:^. (7) 



The diffusion term violates geometric scaling, but is essential to describe the data by mimicking the effect of DGLAP 
evolution. As r decreases, 7off increases from 7s towards the DGLAP value of 1. However, j^g — > 00 in the limit that 
r ^ 0, in disagreement with the expected colour transparency (A/" ^ r"^ as r ^ 0), but this unphysical behaviour has 
been argued to have no influence on the fit since the contribution to cr^ot from very small dipoles is negligible . 

The original CGC fit (and the subsequent fit including charm quarks [l|) fixed the parameters 7s = 0.63 and 
K = 9.9 at the LO BFKL values. The central fits were obtained with Aq fixed at 0.7 and the other three parameters 
((To, xq and A) were fitted to F2 data with XBj < 0.01 and < 45 GeV^. The results from the fit including 

charm quarks [H are shown in the first line of Table HI Statistical and systematic experimental errors are added in 
quadrature. The quark masses in the photon wave functions were taken to be mu^d,s — 0.14 GeV and rric = 1.4 GeV, 
and the dipole cross section was evaluated at a; = XBj for light quarks and x — a;Bj(l + 4to^/Q^) for charm quarks. 
The contribution from beauty quarks was neglected. The saturation scale was found to drop dramatically with the 
introduction of charm; for exainple, the xq parameter was lowered by two orders of magnitude f]\. However, it has 
recently been shown by Soyez that allowing 7s to vary from the LO BFKL value of 0.63 gives an improved fit 
to the F2 data with a 7s slightly larger than 0.7, close to the value extracted from renormalisation-group-improved 
next-to-leading order (NLO) BFKL kernels. Moreover, the saturation scale is not significantly reduced compared to 
the saturation scale obtained with only light quarks. The results of this fit [llj have already been ap plie d to describe 
inclusive diffractive DIS 0| and exclusive diffractive processes using a <-dependent saturation scale [isj . 

In the second line of Table [J we show the results of a fit to F2 data starting from the parameters favoured by Soyez 
PH (keeping k fixed at 9.9). Indeed, the is slightly improved compared to the first fit with a fixed 7s = 0.63. The 
results are not identical to those of Soyez [ll| due to our lack of beauty quark contribution, the fact that we do not 
include HI data, and our more conservative cuts, but none of these differences has much effect on the parameters 
obtained. (As in earlier analyses P, [13, El, we do not fit the HI data to avoid introducing an extra normalisation 
parameter and because the ZEUS data alone suffice.'^ We do not include data with > 45 GeV^ since the CGC 
model does not include the full DGLAP evolution, which becomes more important at large Q^. We do not include 
data with < 0.25 GeV^ because the form of the CGC dipole cross section is motivated by perturbative QCD, 
although including these few data points has little effect on the parameters obtained.) 

In the third line of Table U we show the parameters of a local minimum which is close to the first fit [1], but has a 
slightly smaller 7s = 0.61. Using a hypothesis-testing criterion [2l|, a "good" fit with N degrees of freedom should 
have a of approximately N ± ^/2N . More precisely, we can calculate the p- value of the hypothesis shown in the 
last column of Table IH defined as the probability, under the assumption of a given hypothesis, of obtaining data at 
least as incompatible with the hypothesis as the data actually observed [12] ■ Then, for example, the hypothesis could 
be considered to be excluded at a 90% confidence level if 1 — p > 0.9. The of the first three fits in Table U can all 
be considered to be "good" according to a hypothesis-testing criterion, therefore it is difficult to say that the fit with 
7s = 0.74 is strongly preferred by the data compared to the fit with 7s = 0.61. 

Recent studies of the numerical solution of the BK equation for a fixed impact parameter have shown that the 
effective anomalous dimension T^g is not a constant at the saturation scale, but approaches a limiting value of 
7s ~ 0.44 in the small-x limit [ij].* The value of 7s — 0.44 [l3| seems to be inconsistent with the phenomenological 



^ More precisely, the usual "anomalous dimension" is really 1 — t^w , rather than 7(,fj itself [loll . 

^ However, new preliminary low-Q^ inclusive DIS data from HI [201 will provide tighter constraints on future dipole model fits. 

* In Ref. [12| the LO BK equation with a fixed coupling as = 0.2 was used, which gives a very large value of A ~ 0.9. It was argued that 
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Is 


BcGc/GeV-2 




Xo 


A 


xVd.o.f. 


p- value 


0.63 (fixed, fl]) 


5.5 


0.417 


5.95 X 10^* 


0.159 


211.2/130 = 1.62 


8.7 X 10-*^ 


0.46 


7.5 


0.558 


1.84 X 10"" 


0.119 


118.7/129 = 0.92 


0.73 


0.43 (no sat.) 


7.5 


0.565 


1.34 X 10"^ 


0.109 


124.3/129 = 0.96 


0.60 


0.54 (high-Q^) 


6.5 


0.484 


3.42 X 10"^ 


0.149 


210.9/159 = 1.33 


3.7 X 10"^ 



TABLE II: Parameters of the b-CGC model, ^ and Q, determined from fits to ZEUS F2 data [Ti,[T3 witfi XBj < 0.01 and 
£ [0.25, 45] GeV^. The first line is the fit ,1] with a fixed 7^ = 0.63, while the second line is the new fit where 7^ is allowed 
to go free. The third line is a fit without explicit saturation of the dipole cross section, that is, the form of ^ for rQs < 2 is 
also taken for rQs > 2 (but without the diflusion term in the exponent). The fourth line is a fit which also includes high-Q^ 
DIS data up to < 650 GeV^ 



values found in the CGC model. Fixing 7^ = 0.44 and allowing the other parameters to go free gives the results shown 
in the fourth line of Table HI this cannot be considered a good fit. It is interesting to check if the same conclusions 
are found in the more realistic impact parameter dependent version of the CGC model. 

III. INTRODUCING THE IMPACT PARAMETER DEPENDENCE 

To introduce the impact parameter dependence into the CGC model we modify ([5]) to obtain the "b-CGC" 
model for the dipole cross section P, [2I] : 



^ = 2^(.,.,5) = 2x -^"(^) ^ ^Q^^^ 



d2b 



l_^-Aln^{BrQ^) . , > 2 



(8) 



where the coefficients A and B in the second line are again given by The saturation scale Qs now depends on 
the impact parameter: 



Qs = Qs[x,b) = ( — 



X 



exp 



-1 



2ScGC / 



(9) 



Instead of the normalisation parameter ctq of the CGC model, we now have the parameter BcGCi which is adjusted 
iteratively to give a good description of the t dependence of exclusive J/il) photoproduction; see Sec. IVII We also 
allow the factor Ao to go free. The fit presented in Ref. ,1] fixed 7^ = 0.63, but was unsuccessful in obtaining 
a good fit, as seen in the first line of Table |TT1 For this reason, Ref. [ll| mostly concentrated on the b-Sat model 
description of exclusive processes, that is, using the Glauber-Mueller dipole cross section with DGLAP evolution of 
the gluon density. Nevertheless, the b-CGC fit presented in Ref. [l| has already been applied to calculate electroweak 
deeply virtual Compton scattering 124^ . quarkonium photoproduction in coherent hadron-hadron interactions [2^ . 
the saturation scale in large nuclei [26|, and diffractivc structure functions for both protons and nuclei 27]. 

We find that allowing 7^ to vary in addition to the other parameters dramatically improves the description of the 
F2 data, as seen in the second line of Table HIl Moreover, the optimum value of 7s = 0.46 is close to the expected 
value of 0.44 obtained from the numerical solution of the BK equation [l^. The fact that the optimum value of 7^ 
is quite different in the CGC and b-CGC models, and that the value in the b-CGC model is closer to the theoretical 
expectation, may be attributed to the more realistic modelling of the impact parameter dependence in the b-CGC 
model compared to the factoriscd b dependence implicit in the CGC model. However, the value of A = 0.119 obtained 
from the fit is lower than the perturbatively calculated value of A ~ 0.3 |2^, and suggests that the saturation scale 
comprises significant nonpcrturbative dynamics. 

To examine the importance, or otherwise, of having a unitarised dipole cross section, we performed a fit without 
explicit saturation of the dipole cross section, that is, the form of ^ for rQs < 2 is also taken for rQs > 2 (but 
without the diffusion term in the exponent). The results of this fit are shown in the third line of Table HIl The fit is 
only slightly worse, and the parameters only slightly different, than the main fit presented in the second line of Table 
im This suggests that the HERA data fitted are largely insensitive to the presence of saturation. 



the introduction of a running coupling, giving a smaller value of A, was questionable given that the BK equation is LO in as, and that 
the running coupling did not have a significant effect on the results obtained for the anomalous dimension. 
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Taking the parameters given in the second hne of Table |TT1 and releasing the upper cut on gives a of 381 
for the 163 data points. For comparison, the "b-Sat" model gave a of 193 for the same data [l|. Performing a fit 
improves the to 211 with parameters given in the fourth line of Table El Since the b-CGC model is not justified 
for large Q^, we favour the fit with the more conservative cut of < 45 GeV^, that is, with parameters given in the 
second line of Table [III 

IV. IMPACT PARAMETER DEPENDENT SATURATION SCALE 

It is customary to define a saturation scale Qs, that is, the momentum scale at which the dipole-proton scattering 
amplitude Af becomes sizable such that nonlinear effects start to become important. There is no unique definition of 
Qs and various choices are used in the literature. Following Refs. [l|, ll3| we define the saturation scale Q| = 2/r|, 
where the saturation radius rs is the dipole size where the scattering amplitude 

7V= 1-e-^ -0.4. (10) 

The saturation scale Q| — 2/r| defined by pH)) coincides with Q^(x) = 1/Rl{x) = (xq/x)^ GeV^ in the GBW model 
0, However, for the CGC ([5]) and b-CGC ([8]) models, the saturation scale Qs defined by pO)) differs from the 
parameter Qs- Note that we use uppercase S and lowercase s to distinguish between these two scales. The saturation 
scale Qs is the quantity we shall compute and compare for the different dipole models. 

In Fig. [T] we show the impact parameter dependent saturation scale Q^ for the b-Sat model [ij (solid lines) and 
the present b-CGC model (dotted lines) at 6 = 0,1,2,3 GeV~^. The saturation scale Q^ is strongly dependent on 
the impact parameter b. We also show Q^ for the two candidate solutions for the 6-independent CGC model (dashed 
lines) with 7s — 0.74 and 7s — 0.61. 

In Fig. [2] we show the b dependence in the b-Sat and b-CGC models of the total j*p cross section for representative 
values of x and Q^ of the HERA F2 data included in the fits. The b dependence is approximately factorised from the 
other kinematic variables in the b-Sat model [l|, unlike the case in the b-CGC model. For both models the median 
values of b probed in the total 7*p cross section are in the range 2-3 GeV^"'^. Therefore, from Fig. [1] the solution 
of the 6-independent CGC model with 7s — 0.61 has a similar Q| to the b-CGC model at some average value of b, 
and this is much lower than the Q| of the solution with 7s = 0.74 favoured by Soyez [ll|. Although the b-Sat and 
b-CGC parameterisations of the dipole cross section have very different theoretical motivations, both give a similar 
saturation scale (5| < 0.5 GeV^ for x > 10~^ for the most relevant impact parameters b — 2-3 GeV~^. 

The saturation aspect of the dipole models has been emphasised since the initial GBW investigations 0, Q . How- 
ever, as the dipole models have become more sophisticated, with the introduction of DGLAP evolution, heavy quarks, 
and impact parameter dependence, the saturation scale has lowered and saturation is no longer crucial in describing 
the HERA data. Two-component Regge- motivated dipole models (for example, Ref. [l^) find some preference for 
saturation, but these do not include DGLAP evolution or impact parameter dependence. Indeed, temporary mod- 
ifications of the b-Sat and b-CGC models, so that the dipole cross sections no longer unitarise at large dipole sizes 
and so only the single-Pomeron exchange is present, give comparable fits with only a slight change in the parameters 
required, as seen for the b-Sat model in Ref. 13] and for the b-CGC model in the third line of Table |TT1 Therefore, 
the fact that a model incorporating saturation is successful in describing the data should not be construed as meaning 
that there are large saturation effects present in the data. 

Of course, from the theoretical point-of-view, the dipole model is only self-consistent if ^-matrix unitarity is imposed, 
that is, the scattering amplitude Af cannot take on values greater than one. This feature is necessary to derive relations 
between the inclusive cross section and cross sections for diffractive processes. Moreover, dipole models incorporating 
saturation fitted to HERA data may be extrapolated to very low x (for example, at the LHeC [2^) and to predict 
cross sections for nuclear collisions where the saturation scale is enhanced by A^^^ [26| (for example, at a future 
electron-ion collider (sojV In these situations, multi-Pomeron exchange may become important and extrapolation 
based on single-Pomeron exchange would be unreliable. 

V. DESCRIPTION OF THE LONGITUDINAL AND HEAVY FLAVOUR STRUCTURE FUNCTIONS 

We now compare the predictions of the present b-CGC model, and also those from the b-Sat model [l|, with 
published HERA data on the longitudinal structure function, Fl{x,Q-^), the charm structure function, F|°(a::, Q^), 
and the beauty structure function, i^|''(x, Q^). These predictions can easily be obtained by taking the appropriate 
contributions to the total j*p cross sectionjip]) . 

A first measurement has been reported [31'| of the longitudinal structure function Fl{x, Q^) for Q^ S [12, 90] GeV^ 
based on data taken in the last few months of HERA running, when the proton beam energy was lowered from the 
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FIG. 1: The impact parameter dependent saturation scale Q| = 2/r|, where rs is defined as the solution of (|10|l . at b = 0, 1, 2, 3 
GeV"'^, found in the b-Sat and b-CGC models (second line of Table |TI]). We also show the saturation scale from the two 
candidate solutions of the fo-independent CGC model in Table|I]with 7^ = 0.74 (upper dashed line) and 7s = 0.61 (lower dashed 
hne). We indicate the HERA kinematic limit at y ~ Q^/ixs, s) = 1 with s ~ iEeEp = 90200 GeV^. 



nominal value of 920 GeV to values of 460 GeV and 575 GeV. In Fig. [3] we show these data [3l| compared to predictions 
calculated using ^ with three different dipole cross sections: the b-Sat model l], the present b-CGC model, and the 
GBW fit including charm quarks from Ref. We also show the predictions calculated at next-to-next-to- leading 
order (NNLO) in the conventional coUinear factorisation framework using the MRST 2006 NNLO parton distributions 
[3^ . including the uncertainty band obtained from the 30 alternative eigenvector parton distribution sets. 

In the region of the existing HI data ^31j the three dipole model predictions all agree well with the data and 
also with the NNLO calculation. However, we also show in Fig. [3] the predictions at lower < 10 GeV^. We 
extrapolate the values of x accordingly, by fitting the x values of the HI data as a power law function of Q^, obtaining 
X = (1.09 X 10-5)((52/q2)1-28 ^i^i^ Qo = 1 GeV^. Below 10 GeV^ the dipole model predictions deviate from the 
NNLO predictions and lie well outside the uncertainty band. There is little sensitivity to the detailed form of the 
dipole cross section, as shown by the consistency between the three predictions. The dipole model calculations contain 
some effects from small-x resummation and higher-twist contributions that are not included in the fixed-order coUinear 
factorisation approach, which is known to be perturbatively unstable at low x and Q^. A measurement of at low 
X and would therefore play an important role in discriminatin g be tween the different theoretical approaches. If 
not at HERA, such a measurement could be possible at the LHeC |29| or at a future electron-ion collider [30j . 

In Fig. |3]we show the b-Sat and b-CGC predictions for the charm structure function F2'^{x,Q^) compared to all 
available ZEUS [H, [H, HI] and HI [H, [13, [38,, ^20,] data. Both models give an equally good description of the data, 
even for large x > 0.01 and large where the b-CGC model might be expected to fail. 
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FIG. 2: The 6 dependence of the total 7> cross section, (JtV, for = 0.4, 4 and 40 GeV^ with x = 10"*, 10"^ and 10" 
respectively, in the b-Sat model (left) and the b-CGC model (right). The median values of b are all between 2 and 3 GeV"^. 




FIG. 3: Predictions for the longitudinal structure function, Fl(x, Q^), using three different dipole model cross sections compared 
to HI data 31]. We also sh ow p redictions calculated at NNLO in the coUinear factorisation framework using the MRST 2006 
NNLO parton distributions 




FIG. 4: Predictions for the charm structure function, F2'^(x,Q^), using two different dipole model cross sections compared to 
data from ZEUS [3^. 13^. [SSj and lil 36, 37, 38, 39]. The bins used correspond to the data points in Ref. j34| ; the other 
data points have been shifted to these values using the b-Sat model predictions. 
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FIG. 5: Predictions for the beauty structure function, F2' {x^Q^), using two different dipole model cross sections compared to 
data from HI 
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b-Sat model 
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FIG. 6: The 7*p scattering amplitude in the b-Sat model (left) and the b-CGC model (right). For inclusive DIS, we have 
-E = 7*, a; = a;' <C 1 and p = p' . For exclusive diffractive processes, such as vector meson production [E = V = J/^}j,4',p) or 
DVCS {E — 7), we have x' <^ x ^ 1 and t = {p — p')^. The dipole cross section is unitarised assuming eikonalisation of the 
6-dependent DGLAP-evolved gluon density in the b-Sat model (left), or via the fan diagrams resummed by the BK equation 
in the b-CGC model (right). 



In Fig. [5] we show the predictions for the beauty structure function F|''(a;, Q^) using a bottom quark mass of 
TOf, = 4.5 GeV. (There is some sensitivity to this value, especially at low Q^.) Although the available data have large 
statistical uncertainties, the difference between the b-Sat and b-CGC model predictions increases as increases, 
such that the b-Sat model predictions are favoured over those from the b-CGC model in the higher bins. 



VI. DESCRIPTION OF EXCLUSIVE DIFFRACTIVE PROCESSES AT HERA 



We now check that the updated b-CGC model, with the parameters give n in the second line of Table HTl describes 
the main features of exclusive diffractive vector meson (J/-!/! fid, l4ll. |42||. (f> [i^ and p [H, H^) production and 
deeply virtual Compton scattering (DVCS) [H, \^ [H, [49,] at HERA, and compare the predictions with the b-Sat 
model [l|. The main concepts of the b-Sat and b-CGC models are illustrated in Fig, p The new precise data on p 
electroproduction from ZEUS [ZF] and new data on DVCS from HI [48*1 and ZEUS [ij which were not available at 
the time of Ref. [1], provide tighter constraints on the models. The data have been measured as a function of the 
photon virtuality, Q^, the 7*p centre-of-mass energy, W, and the four-momentum transfer at the proton vertex, t. In 
addition, the ratio of the cross sections for longitudinally and transversely polarised photons has been measured. It is 
therefore a significant challenge for an essentially parameter-free model to describe all features of the available data. 

The formalism for calculation of exclusive diffractive processes from the 6-dependent dipole cross section was 
discussed in detail in Ref. [ij. Here we give only the final formulae. The differential cross section for the exclusive 
process 7*p Ep, where E — V,j, is [l] 



1 

167r 



A 



7*p — >Ep 
T.L 



(11) 



where the scattering amplitude 



47r 



d2b 



(12) 



For DVCS, the amplitude involves a sum over quark flavours f — u^d, s, c. The ratio of the real to imaginary parts 
of the scattering amplitude, /3, is calculated using 



13 — tan 



ttA 



with A 



a In (y^j'i:^^^ 

d\n{l/x) 



(13) 



The factor in (fTTj) accounts for the skewedness effect, that is, x 7^ x' in Fig. [6l and is calculated using [SC 



(14) 
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10 10^ 1 10 1 10 10^ 

+ Mj,^ (GeV^) q2 + (GeV^) + (GeV^*) 

FIG. 7: The total cross section cr vs. {Q^ + My) for exclusive J/^jj [iol. l4ll. . [i^l and p [45l | meson production compared 
to predictions from the b-Sat and b-CGC models using the "boosted Gaussian" vector meson wave function. The ZEUS J/tp 
photoproduction point is taken from Table 1 of Ref. [4C1( ]. from the muon decay channel with W = 90-110 GeV. 

with A calculated as in (jl3p .^ The forward overlap functions between the photon and exclusive final state wave 
functions in (fT^ . are given in Ref. [l]]. We use the "boosted Gaussian" vector meson wave functions 

[Sllls^lssj. which were found to give the best description of data in Ref. [ij. We take x — xbj (1 + My/Q^) for vector 
meson production. For DVCS we take x = xbj for the light quark contributions and x — XBj(l + 4to^/(3^) for the 
charm quark contribution. 

In Fig. [7] we show the description of the + My dependence of the vector meson data. In Fig. [5] we show the W 
dependence, and in Fig. [9] we show the dependence of 5, where a cx W^. In Fig. [10] we show the dependence of 
Bd, where dtr/dt cx ex\i(— B D\t\) ■ In Fig.[TT]we show R = OhjaT vs. Q^, W and t, for p meson electroproduction. In 
Fig.[T2]we show the effective Pomeron trajectory ap(i) vs. |f|, where ap(t) is determined by fitting dti/di cx I^^["p(*)^il . 
In Fig.[Tn]we show the i-slope parameter Bn vs. W . Finally, in Figs . [HI and [T51 we show the Q^, W and t dependences 
for DVCS. 

In general, both the b-Sat and b-CGC dipole models, using the "boosted Gaussian" vector meson wave function, 
describe almost every feature of the available data. The b-CGC model gives a better description of ap, where 
Q!p(i) = ap(0) -I- apt, as noted in Ref. [1|], as seen from the slope of the predictions in Figs. [T2] and [T3l One major 
problem, again as noted in Ref. [l[, is with the ratio R = gl/gt for p meson electroproduction at large Q^; see 
Fig. [TT] The flat behaviour of i? as a function of |t| seen in the data implies that the t-slope parameter Be, is the 
same for both transverse and longitudinal photon polarisations. This surprising behaviour contradicts both the model 
predictions and the general intuitive picture, where the transversely polarised cross section is dominated by larger 
dipole sizes than the longitudinally polarised cross section, leading to Bd.t > Bd.l- A detailed examination of the 
transversely polarised p meson wave function is needed to resolve this issue. 

The W dependence of J/^p photoproduction is much better described by the b-Sat model than by the b-CGC 
model. (The description with the b-CGC parameters given in the fourth line of Table [III where high-Q^ F2 data were 
also included, is not any better.) Indeed, the W dependence, or the power S, is a test of the evolution of the dipole 
cross section (or generalised gluon density) and provides a powerful discriminator between the different models. In 
Fig. I16f a) we show the r dependence of the imaginary part of the scattering amplitude for J/ip photoproduction. 
The distributions are peaked at r « 1.1-1.3 GeV~^, cf. the "scanning radius" [5l| of r w = 1.9 GeV~^. 

Going from W = 50 GeV {x = 4 x IQ-^) to = 300 GeV {x = IQ-*), the b-CGC model amplitude goes from 
being larger than the b-Sat model amplitude around the peak to being smaller. This behaviour can be traced to 



^ In Ref. T] the factor Rg in JTTJ was omitted for the b-CGC model. For the b-Sat model, the factor Rg instead multipHes the gluon 
density inside the dipole cross section. 
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FIG. 8: The total cross section a vs. W for exclusive J/ip [iol. l4ll. |43| . (f> [43| and p p4. l4a| meson production compared to 
predictions from the b-Sat and b-CGC models using the "boosted Gaussian" vector meson wave function. The ZEUS J/ip 
data points [i^, |4j] have been scaled to the HI values using the dependence measured by ZEUS of the form 



[41l |. The HI p data points I4J| have been scaled to the ZEUS Q values [45j using the Q dependence 



measured by HI of the form a oc (Q -|- My) 
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FIG. 9: The power S vs. (Q^ + My), where S is defined by fitting a oc W\ for exclusive J/V' [13, Sll, SI, <^ [H and p [H 
meson production compared to predictions from the b-Sat and b-CGC models using the "boosted Gaussian" vector meson wave 
function. 



the 6-integrated dipole cross section in the relevant r region, shown in Fig. [TBf b) at two different values of x. For 
comparison, we also show the GBW fit, including charm quarks, from Ref. pj. To examine the differences between 
the three parameterisations of the dipole cross section in more detail, we can also compute the effective anomalous 
dimension, 7cff = 91ntTgg/91nr^, shown in Fig. [TBTc). For small dipole sizes, the value of 7eff for the b-Sat (GBW) 
dipole cross section approaches a limiting value of approximately (exactly) 1. On the other hand, the value of 7cff for 
the b-CGC model diverges for small r and large x\ see ([7]). 
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FIG. 10: The t-slope parameter Bd vs. {Q^ + My), where Bn is defined by fitting da/dt oc exp{—BD\t\), for exclusive J/ip 
[40l . I41I . [i^ . (f) [i^l and p [il . fisj meson production compared to predictions from the b-Sat and b-CGC models using the 
"boosted Gaussian" vector meson wave function. 




Although there are relatively small differences between the b-Sat and b-CGC dipole cross sections shown in Fig. fTBf b) 
around r ^ 1 GeV~^, the J/ip data can discriminate between them, as seen more clearly in the description of the 
6 parameter for J/iJj photoproduction shown in Fig. [HI Note that the skewedness and real part corrections in (jlip . 
not included in Fig. [TBI enhance the relative difference between the b-Sat and b-CGC dipole cross sections by about 
a factor 1.2 due to the larger values of A obtained in the b-Sat model. A similar difference between the b-Sat and 
b-CGC model predictions is found in the values of 6 for DVCS shown in Fig. [TU however, the HERA data on DVCS 
are not yet precise enough to distinguish between the two models. 

We conclude that the b-CGC model, and other similar models derived from (O, work well for the designed purpose 
of extracting the saturation scale and describing low-to-moderate DIS, and also for describing exclusive diffractive 
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The total cross section a vs. and W, and the power 5 vs. Q^, where S is defined by fitting a oc , for DVCS 
I47I. l48l. I49II . compared to predictions from the b-Sat and b-CGC models. 



p and 4> meson production and DVCS. However, such models fail for observables sensitive to relatively small dipole 
sizes, such as J/V' photoproduction and high-Q^ DIS, where the b-Sat model [H with explicit DGLAP evolution 
performs better. Note that the failure of the b-CGC model to describe J/ip photoproduction does not indicate that 
the model should also fail for inclusive DIS at ~ ^-^J/jp' si^*^^ exclusive diffractive processes are known to be 
dominated by smaller dipole sizes than inclusive DIS at the same Q^; see, for example, Ref. [l^. 



VII. CONCLUSIONS 

We have extended the CGC dipole model of lancu, Itakura and Munier [l3| to include the impact parameter 
dependence. The value of the anomalous dimension at the saturation scale, 7s = 0.46, determined from fitting to the 
total 7*p cross section at HERA, is close to the value of 0.44 obtained from numerical solution of the BK equation 
p^2j . However, the value of A = 0.119 is lower than might be expected from a perturbative calculation [28j. The 
impact parameter dependent saturation scale, Ql, is generally < 0.5 GeV^ in the HERA kinematic regime for the 
most relevant b ~ 2-3 GeV""'^. This is in agreement with previous findings from the b-Sat model P, that is, the 
Glauber-Mueller dipole cross section with DGLAP evolution of the gluon density. We have shown that the data do 
not show a strong preference for the solution of the CGC model presented by Soyez [ll| , which has a large saturation 
scale Ql ~ 1 GeV^ at a; = 10~^. It would be interesting to check whether the present b-CGC dipole cross section is 
compatible with data from RHIC. The dipole model predictions have been successfully compared with the first direct 
measurement of the longitudinal structure function from HERA [sij , and also with existing HERA data on the charm 
and beauty structure functions. 

Although a leading-twist collinear factorisation theorem has been proven _54i| by which the amplitude for exclusive 
meson production can be expressed as a convolution of generalised parton distributions (GPDs) with hard-scattering 
kernels and meson distribution amplitudes, there are a number of practical problems with this approach (more 
progress has been made for DVCS (55|). Firstly, the GPDs are not well-known and one needs to rely on models where, 
for example, the GPDs are written in terms of the usual parton distributions of the proton. (However, the same 
criticism could be applied to the modelling of the dipole cross section.) Secondly, the NLO corrections have been 
calculated and found to be huge at small x [56| , implying that small- a; resummation will be needed to achieve a stable 
result. Therefore, the colour dipole approach discussed in this paper, or the related /cj -factorisation approach (see, for 
example, Refs. [13, HI]), provide a complementary way to describe exclusive diffractive processes. These approaches 
include some effects from small- a; resummation, and also corrections due to the transverse momentum of the partons 
entering the hard-scattering subprocess, which are neglected in the collinear factorisation approach but are known to 
be substantial [1^. Of course, the weakness of the dipole picture with respect to the collinear factorisation approach 
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FIG. 16: (a) The r dependence of the amphtude for exclusive J/t/j photoproduction at t = (temporarily omitting the 
skewedness and real part corrections), (b) the 6- integrated dipole cross sections divided by r^, and (c) the effective anomalous 
dimension, jcS = dlnaqq/dlnr'^ . 



is that it is not currently known how to systematically improve the dipole picture beyond LO. 

The b-CGC model presented in this paper has been shown to provide a good description of exclusive diffractive p 
and meson electroproduction and DVCS, but performs less well for observables sensitive to relatively small dipole 
sizes, such as exclusive diffractive J/^p photoproduction and the beauty structure function at high Q^. In these cases, 
the b-Sat model [l| with explicit DGLAP evolution performs better. The W dependence of exclusive diffractive 
processes provides an important discriminator between different dipole model cross sections. It would therefore be 
interesting to measure exclusive diffractive processes with greater precision than at HERA, as may be achievable at 
the LHeC [1^ or with a future electron-ion collider [soj . 

It should be borne in mind that the dipole picture is not exact and that there are a number of assumptions and 
approximations made in its formulation [sl,!!^]. Nevertheless, despite the de-emphasis of the saturation aspect of the 
dipole models, the framework provides an intuitive and economical description of a wide variety of small-x processes. 
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